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Neutral-meson mixing is loop suppressed in the Standard Model, lead¬ 
ing to the possibility of enhanced sensitivity to new physics. The uncer¬ 
tainty in Standard Model predictions for S-meson oscillation frequencies is 
dominated by theoretical uncertainties within the short-distance S-meson 
hadronic matrix elements, motivating the need for improved precision. In 
D-meson mixing, the Standard Model short-distance contributions are 
further suppressed by the GIM mechanism allowing for the possibility of 
large new physics enhancements. A first-principle determination of the 
D-meson short-distance hadronic matrix elements will allow for model- 
discrimination between the new physics theories. I review recently pub¬ 
lished and ongoing lattice calculations of hadronic matrix elements in B 
and T)-meson mixing with emphasis on the Fermilab lattice and MILC 
collaboration effort on the determination of the B and T)-meson mixing 
hadronic matrix elements using the methods of lattice QCD. 
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1 Motivation 


In the Standard Model (SM), neutral-meson mixing is loop suppressed due to the 
absence of a tree-level flavor changing neutral current, resulting in greater sensitivity 
to physics beyond the SM. Historically the neutral Kaon and S-meson system had 
inferred the charm quark and top quark masses respectively before direct production. 
In the pursuit of indirect detection through precision physics, current and planned 
experiments from LHCb, BELLE II and BES III promise continual effort in the flavor 
frontier. Therefore, in tandem with current experimental measurements, and in antic¬ 
ipation of future results, we calculate the hadronic matrix elements that participate 
in neutral B and D-meson mixing using the methods of lattice QCD. 

2 Neutral-meson mixing phenomenology 
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(a) Short-distance diagram. The external 
state H indicates a general heavy-light 
neutral-meson. In the SM, the internal 
quarks q are up-type for B- and down- 
type for T)-meson mixing. 
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(b) Long-distance diagram. The exter¬ 
nal state H indicates a general heavy-light 
neutral-meson. Dominantly a SM process, 
this diagram proceeds via numerous multi¬ 
hadron and resonant intermediate states. 


Figure 1: Contributions to neutral-meson mixing. 


Mixing of the neutral-meson states in general receives contributions from short- 
and long-distance diagrams. The Feynman diagram of the short-distance contribution 
is sketched in Fig. The initial and final state hadrons are described by strong 


dynamics at the MeV scale of Aqcd; in the SM, mixing occurs at the weak scale and 
is in comparison, a short-distance interaction. For the S-meson, the SM box-diagram 
is enhanced by the top-quark by a factor of ^ 20, and is the dominant 

process describing S-meson mixing. Diagrams mediated by the up quark and charm 
quark are subdominant and are further suppressed by the GIM mechanism. For 
D-meson mixing, the diagram mediated by the bottom quark is CKM suppressed 
by r\j 0.2^°, while diagrams proceeding via the down quark and strange quark 
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are GIM suppressed. Therefore, D-meson mixing in the SM receives subdominant 
short-distance contributions. 

Neutral-meson mixing receives long-distance SM contributions mediated by inter¬ 
mediate hadronic states, as shown in Fig. [TO For 5-meson mixing, the top-quark 


can not hadronize, while hadronization via the up quark and charm quark are CKM 
suppressed by at least ~ 0.2^, and therefore is a subdominant contribution. On the 
D-meson side, while the bottom quark receives the same CKM suppression as seen 
in the short-distance diagram, intermediate hadrons with down quarks and strange 
quarks are only doubly Cabbibo suppressed. Several model dependent estimates 
have shown that these long-distance contributions may reach the observed level of 
D-meson mixing [1]. However, the exclusive determination of the long-distance dia¬ 
gram is complicated by many-body multi-channel and resonant intermediate states, 
while inclusive calculations rely on the poor expansion parameter of Aqcd/?^c- In 
both cases, estimates are beset by large hadronic uncertainties. 

Beyond the SM, neutral-meson mixing proceeds via short-distance diagrams since 
the scale of new physics is expected to be above the electroweak scale. New physics 
may facilitate 5-meson mixing, and can be uncovered by observing a discrepancy 
between observed 5-meson mixing and the SM prediction. 5-meson mixing has am¬ 
ple room to proceed dominantly through physics beyond the SM due to a weakly 
constrained SM description. Since the phase between the SM and new physics contri¬ 
bution is unknown, one assumes that the SM contribution is subdominant and derive 
bounds on new physics [2] . Therefore, the lattice community is motivated to provide 
determinations of all hadronic matrix elements pertinent to 5- and 5-meson mixing 
at the level of precision comparable to experimental observation. 

All hadronic matrix elements involved in neutral-meson mixing can be expressed 
in a basis of five 4-quark operators, 

Oi = O2 = O3 = W" , ( 1 ) 


O4 = 


O5 = 


( 2 ) 


where 4' and ip denote the heavy and light valence quarks of the 5 or 5-meson, L and 
R are left and right-handed projection operators, and the Greek indices label color 
degrees of freedom. Operator Oi has the SM charged current structure, and is directly 
related to the matrix element of the box diagram. Operators O 2 and O 3 couple only to 
left-handed quarks, and enter in the calculation of the width difference [3]. Operators 
O 4 and O 5 couple to right-handed quarks and are in general needed for new physics 
models. Therefore, on the lattice we need only calculate five matrix elements for each 
neutral-meson to describe all hadronic interactions. 
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(a) MILC asqtad gauge ensembles. Colors la¬ 
bel ensembles with different lattice spacing, 
and darker hues correspond to lighter pion 
masses. The cyan star marks the physical 
pion mass in the continuum limit. Ensem¬ 
bles include 592 to 2258 uncorrelated Monte 
Carlo samples indicated by the size of the 
markers. 


(b) Light-quark propagators. Statistically 
correlated light-quark propagators of differ¬ 
ent valence masses are obtained at each 
gauge configuration. The black diagonal line 
passes through mgea = '^valence- A total of 
102 data points are available per matrix ele¬ 
ment to guide the chiral-continuum extrapo¬ 
lation. 


Figure 2: Dataset used for the analysis. 


3 Lattice calculation 

The Fermilab lattice collaboration calculates hadronic matrix elements for Bd, Bg, 
and D using a large set of the MILC collaboration’s asqtad gauge ensembles [1]. The 
light valence quark propagators are generated using the asqtad action, and the heavy- 
valence quark propagators use the Sheikoleslami-Wohlert action with the Fermilab 
interpretation [H] . A detailed list of ensembles used is tabulated in Ref. [H] . Illustrated 
in Fig. 1^ the wealth of data over a wide range of quark masses and lattice spacings 
allows us to control and estimate all systematic errors. A summary of the inclusion 
of various systematic errors is detailed in Ref. [3 IE]. We perform a combined chiral- 
continuum extrapolation over the entire data set and obtain the matrix elements at 
the physical point ng, as illustrated in Fig. Preliminary results reveal that the 
Fermilab lattice calculation for matrix elements have errors ranging from 5% to 15%; 
a detailed error breakdown is arranged in Tab. 
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valence light quark mass [ri] 


Figure 3: Chiral-continuum extrapolation. Extrapolation to the physical point for 
the D-meson mixing matrix element of operator Oi. The data points share the 
color coding of Fig. to reflect the gauge ensemble in which they are evaluated on. 
The continuum extrapolation is marked by the cyan band, the yellow line marks the 
physical up-quark mass. 



1 {HmH) 

(771021^-) 


(771041^-) 

{HmH) 1 

Bd 

7.7% 

8.4% 

15.2% 

6.7% 

8.9% 

Bs 

5.5% 

5.8% 

10.8% 

5.1% 

6.6% 

D 

8.7% 

4.7% 

6.1% 

5.0% 

9.1% 


Table 1: Preliminary relative error for the B^, Bg, and D-meson hadronic matrix 
elements. 
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4 Phenomenological implications 

4.1 B-meson mass splitting 

In the SM, the S-meson mass splitting is given by, 

where the factor in the parenthesis and riB are know electroweak contributions with 
QCD corrections. In combination with (^Bq\Oi\Bq'j, the theoretical prediction can 
then be compared with the experimental observation of AMq. 

The matrix elements are often expressed in terms of the bag parameter, 

{H\Oi\H)=CMfHB^H\ (4) 

where i labels one of five 4-quark operators, Mh is the mass of the heavy-light meson 
H, fn is the meson decay constant, B^^ is the bag parameter, and Ci is a factor that 
normalizes B^^ in the vacuum saturation approximation. 

The latest experimental observations of the S-meson mass differences are 

AMd = 0.510(3) ps-^ AMs = 17.761(22) ps-\ (5) 

with respective errors of 0.6% and 0.1% [H]. The current FLAG-averaged lattice 
determination of the bag parameters Bb^^ comes entirely from HPQCD with 8% and 
5% errors [laiiH]. The Fermilab lattice collaboration will provide an independent 
determination of the bag parameters, however preliminary results suggest similar 
errors. There are ongoing efforts to reduce the hadronic uncertainty with hints of 
sub-percent precision for ratio quantities [T8] . 

4.2 Unitarity triangle constraint 

The S-meson mass splittings constrain one side of the unitarity triangle through the 
following relationship 

bid _ ^ / AMdMss 

where is the ratio of the Bg and Bd matrix elements for operator Oi. On the lattice, 
^ is better determined than the bag parameters since taking the ratio partially cancels 
statistical and systematic errors. The current FLAG average of ^ comes entirely 
from Fermilab lattice with 5% error da E]. Preliminary results from Fermilab 
lattice reduce this error down to 1%, largely as a result of enlarging the data set and 
controlling the systematic error stemming from wrong-spin taste-mixing effects [H]. 
The projected improvement to the constraint on the unitarity triangle is illustrated 
in Fig. m where the pink bounds shrink to the red bounds. 
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Figure 4: Unitarity triangle constraints from semileptonic B decays and S-meson 
mixing at two standard deviations. Pink curves are evaluated at the FLAG average ^ 
value, red curves are evaluated at the FLAG average central value with the Fermilab 
lattice preliminary error. All other parameters are taken from Ref. m- 

4.3 D-meson mass splitting 

Assuming that D-meson mixing proceeds dominantly through BSM processes, the 
mass splitting is given by, 

AM = Ci,p {D (7) 

where CJjp is the model-dependent short-distance coefficient. Supplied with exper¬ 
imental observation for AM, and lattice calculations for (Oi), parameters in new 
physics models can be constrained. One may also provide model independent con¬ 
straints to new physics following Ref. [E] and [Tf)] . 

The current experimental measurement for the T)-meson mass splitting is 

AM = 0.0044(20) ps"^ (8) 

with a 45% error [TT] . Gurrently there is one published result on the D-meson mixing 
hadronic matrix elements from the ETMG with 2 dynamical sea quarks [16] . Results 
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from Fermilab lattice will be the first 2+1 dynamical sea quark calculation and pre¬ 
liminary results, shown in Fig. suggest the calculations are statistically consistent. 
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Figure 5: D-meson mixing matrix elements from ETMC and preliminary Fermilab 
lattice results. Bag parameters from ETMC [Hj are converted to matrix elements 
with the D-meson decay constant from Ref. m. 


Planned experiments for the next decade indicate that the error is expected to be 
reduced to the 10% level, on par with lattice uncertainties of Tab. [T] 


5 Conclusions and outlook 

We are calculating Bg, and D-meson hadronic matrix elements with 5% to 10% 
relative error, and the ratio quantity ^ at the 1% level. For the B system, hadronic 
quantities are currently the dominant source of uncertainty and ongoing efforts from 
the lattice community aim to obtain errors commensurate with experiment. Un¬ 
certainty on the H-meson hadronic matrix elements are comparable to future experi¬ 
mental observation over the next decade. However, uncontrolled hadronic uncertainty 
exists in the SM long-distance contribution. 
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